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ABSTRACT 


The  theoretical  acoustic  behavior  of  an  underwater  bubble 
screen  was  examined  using  Rayleigh  reflection  theory.   A 
microcomputer  model  simulates  the  acoustic  impedance  mismatch 
at  the  interfaces  of  an  ideal,  bubble  screen.   A  sensitivity 
analysis  indicates  that  the  angle  of  incidence  of  sound 
energy  and  the  speed  of  sound  in  the  layer  are  the  most 
important  screen  properties  for  predicting  the  insulating 
capability  of  a  bubble  screen.   In  the  neighborhood  of 
frequencies  for  which  the  screen  thickness  is  an  integral 
number  of  half  wavelengths,  the  interference  results  in  a 
reduced  reflection  coefficient  and  a  corresponding  increase 
in  transmission  through  the  screen.   So  that  for  a  broad 
band  spectrum  wide  enough  to  cover  a  number  of  such 
frequencies,  the  attenuation  to  be  expected  exceeds  10  dB 
only  over  about  90  percent  of  the  spectrum.   The  interest 
for  this  work  came  from  the  need  for  a  noise  insulating 
screen  at  the  Carr  Inlet  Acoustic  Range. 
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I o    INTRODUCTION 

Acoustic  trials  of  ships  and  submarines  are  conducted  by 
Puget  Sound  Naval  Shipyard  (PSNS)  at  the  Carr  Inlet  Acoustic 
Range  in  the  state  of  Washington.   According  to  Mr o  John 
Kriebel  of  PSNS  (personal  communication)  the  measurements 
are  often  contaminated  by  noise  generated  by  ship  traffic 
in  adjacent  waters.   One  method  which  may  offer  noise  isola- 
tion is  a  bubble  screen  located  between  the  range  and  the 
offending  noise  source.   Mr ^  Kriebel  suggested  that  a 
feasibility  study  be  conducted. 

This  thesis  is  one  part  of  the  study  and  has  the 
objective  of  examining  the  theoretical  benefits  of  a  bubble 
screen,  using  a  computer  model  based  on  Rayleigh  reflection 
theory  for  an  homogeneous  layer.   Recommendations  and  con- 
clusions are  made  based  on  a  literature  search  and  a 
sensitivity  analysis  of  the  computer  modelo   A  companion 
thesis  by  LT  Kelley  examines  the  noise  generated  by  the 
bubbles  themselves. 

The  writer  wishes  to  acknowledge  the  interest  and 
support  from  the  PSNS  Carr  Inlet  Acoustic  Range  during 
this  study. 
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II o    BACKGROUND 

Ao    RESONANCE  BUBBLE  THEORY 

Experimental  studies  of  the  acoustic  behavior  of  bubbles 
in  water  have  been  conducted  for  many  years o   Normally,  when 
a  single  bubble  is  being  studied,  the  resonance  theory  of  gas 
bubble  pulsations  proposed  by  Minnaert  in  1933  [Ref.  1]  is 
considered.   This  theory  shows  that  a  gas  bubble  in  water  is 
capable  of  vibrations  which  become  large  at  an  exciting 
frequency,  f,  in  hertz  given  by 


^•"s  2Kr   V  9 


where , 


1^  =  ratio  of  the  specific  heats  of  the  gas 

P  =  ambient  hydrostatic  pressure  in  pascals 

O  =  density  of  water  medium  in  kilograms  per  cubic 
'  meter 

r  =  radius  of  the  bubbles  in  meters 

Table  I  shows  the  radius  of  a  resonant  bubble  as  a 

function  of  frequency  f  and  depth  z  of  concern  at  Carr  Inlet 

TABLE  I 
Resonant  Bubble  Sizes 


0  m 


50  m 


100  m 


150  m 


200  m 


0,32    m            0.79    m            1,08    m  1.3    m            1,5    m 

0»032m         0.079m         0,108m  0,13m          0.15m 

3  .  2    mm            7  ,  9    mm          1 0  ,  8    mm  1 3    mm            1 5    mm 

0.32    mm         0,79    mm          1,08    mm  1 , 3    mm          1 . 5    mm 

3      3 
where   ^  =   1.03  x  10   kg/ra 

P  =   10   (1  +  0,lz)  pascal 


If  =   1.40 
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Carstenson  and  Foldy  [Refo  2J  examined  certain  acoustic 
properties  of  bubble  screens  and  demonstrated  reasonable 
experimental  agreement  with  the  theory  of  resonance  of 
bubbles  in  water. 

Although  resonance  theory  predicts  very  large  attenua- 
tions of  incident  sound  energy  for  bubbles  in  water,  this 
attenuation  is  a  very  strong  function  of  the  bubble  sizeo 
"75  percent  of  the  energy  lost  by  the  sound  wave  is  due  to 
bubbles  of  radii  within  10  percent  of  the  resonant  radii," 
[Refo  3:  pol3j   A  parameter   used  to  describe  bubble  behavior 
is  its  extinction  cross  section  (^)  /  the  ratio  of  the 
energy  lost  from  the  sound  wave  to  the  intensity  incident  on 
the  bubble.   [Ref.  4j   ^   is  made  up  of  both  absorption  cross 
section  {(JZ  )     and  scattering  cross  section  {(JZ  )     where 

•  ^  =  Q^f  ^ 
For  a  plane  wave  propagating  in  a  bubbly  mixture,  the 
attenuation,  a,  for  bubbles  of  uniform  radius  r  is 
{Ref.  5:  P.615J 

a    =  ^.3^'g^-n(r)        iV  dB  per  /vieter 

where  n(r)  =  number  of  bubbles  per  cubic  meter  of 
radius  r 

Medwin  [Ref,  5j  has  indicated  that  ".^^at  a  resonance,  the 

scattering  and  absorption  cross  sections  of  a  typical  bubble 

at  sea  are  of  the  order  1000  times  its  geometric  cross 

section, " 
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Several  points  should  be  made  here. 

1.  The  sensitivity  of  attenuation  on  bubble  sizes 
requires  a  precise  control  of  bubble  size. 

2.  As  will  be  discussed  later,  even  if  bubbles  of 
radius  greater  than  1  to  2  cm  could  be  produced,  these 
bubbles  will  rapidly  become  unstable  and  break  up. 

3.  A  bubble  will  expand  mono tonically  as  it  rises  to 
accomodate  the  decrease  in  hydrostatic  pressure.   A  bubble 
of  1  mm  radius  at  100  m  will  expand  to  3.2  mm  at  the  surface 
and  its  resonance  frequency  will  vary  from  10  kHz  to  1.0 

kHz  during  its  ascent. 

4.  MacPherson  [Ref.  7J  indicates  that  it  is  difficult 
to  produce  bubbles  of  a  given  size.   Consequently,  it  is 
difficult  to  predict  the  attenuation  obtainable  from  given 
bubble  production  mechanisms. 

B.    RAYLEIGH  REFLECTION  THEORY 

For  the  above  reasons  and  because  of  the  difficulty  of 
producing  bubble  resonance  at  low  frequencies,  this  paper 
proposes  to  view  the  bubble  screen  as  an  homogeneous  layer 
separating  semi-infinite  layers  of  bubble  free  water. 
Rayleigh  reflection  theory  predicts  that  attenuation  of 
sound  energy  by  the  layer  of  bubbles  will  be  a  function  of 
layer  width,  angle  of  incidence,  frequency,  and  the  acoustic 
properties  of  the  bubble  layer.   A  layer  of  constant  width 
will  be  used  to  simulate  the  bubble  screen. 
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A  sensitivity  analysis  will  then  indicate  which 
parameters  are  of  greatest  concern  in  bubble  layer  design. 
The  sound  power  transmission  coe  f  f  icient ,  c?<-t- ,  will  be  the 
measure  used  to  analyze  these  parameters.   cKlis  related  to 
the  sound  power  reflection  coefficient  of  by  the  equation, 

K 

C.   MANUFACTURING  BUBBLES 

In  the  past,  researchers  used  a  variety  of  techniques  to 
control  the  size  of  bubbles  in  screens.   The  methods  attempted 
include  electrolysis,  JRef.  8j  microdisper sers ,  [Ref.  9j 
porous  porcelain  filters,  [Ref.  loj  and  perforated  pipes. 
Exact  bubble  size  was  important  to  determine  acoustic 
characteristics  of  scattering  cross  section,  extinction  cross 
section,  resonance  frequency,  and  other  parameters.   Several 
papers  concerning  bubble  generation  [Refs.  11,  12,  13,  14] 
indicate  that  some  of  the  many  factors  controlling  bubble 
size  for  perforated  pipes  include: 

1 .  orifice  size 

2.  pressure  drop  across  the  orifice 

3.  wetting  of  the  orifice 

4.  contaminants  in  the  media 

5.  angle  of  inclination  of  the  orifice 
5.   rate  of  bubble  formation 

Uncertainties  in  these  factors  are  the  primary  factors 
which  will  prevent  accurate  bubble  size  prediction  for  a 
bubble  screen  installation  on  a  large  scale  such  as  at  Carr  Inlet 

15 


D.    BUBBLE  HYDRODYNAMICS 

Air  bubbles  in  water  exhibit  complex  hydr odynamic 
characteristics.   However,  their  behavior  can  be  categorized 
generally  with  respect  to  the  bubble  size.   Initial  bubble 
motion  has  been  analyzed  by  Walters  and  Davidson 
[Ref s .  15  and  16j  and  Bachhuber  and  Sanford  jRef.  17]. 
Gas  diffusion  from  bubble  affects  bubble  behavior  signif- 
icantly when  the  bubble  diameter  is  small.   The  gas  diffusion 
effect  depends  on  the  rate  of  ascent,  the  radius  of  the 
bubble  and  the  difference  between  the  gas  concentrations 
across  the  bubble  liquid  interface.    [Ref.  18:  p.  714j 
There  are  two  critical  radii: 

R    =  the  radius  above  which  the  bubble  qrows  monotonical ly 
in  size  as  it  ascends;  gas  leakage  is  not  predominant 

R    =  the  radius  below  which  the  bubbles  shrink  and 

collapse  bubbles  whose  radii  lie  between  Ra  and 
Rb  are  unpredictable. 

Two  bubbles  of  the  same  gas  and  size  produced  at  the  same 

depth  can  be  affected  differently  by  surfactants,  causing 

one  of  the  bubbles  to  grow  and  the  other  to  shrink  and 

collapse. 

Bubbles  of  interest  for  a  screen  at  Carr  Inlet  are  in 
the  1-10  mm  diameter  range.   This  size  has  proven  to  be  the 
easiest  to  produce  via  a  perforated  pipe.   This  method  is 
also  by  far  the  most  convenient  and  cost  effective. 

"Hydrostatic  head  has  very  little  influence  on  the 
relation  between  the  rate  of  formation  and  the  size  of  the 
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bubble."   [Ref.  19:  p. 17]   VanKrevelen  and  Hoftijzer 

[Ref.  20:  p.30j  showed  that  "...the  diameter  of  the  bubbles 

is  independent  of  flow  rate  and  increases  with  the  cube  root 

of  the  orifice  diameter."   This  statement  breaks  down  at 

higher  flow  rates  when  gas  jetting  effects  begin. 

Crump  [Ref.  2l]  suggests  a  differential  pressure  of  2-3 

psi  be  established  between  pipe  line  pressure  and  hydrostatic 

water  pressure.   This  should  make  the  desired  bubbles  for 

forming  the  bubble  screen.   A  bubble  of  radius  Imra  will 

expand  monotonically  due  to  decreased  hydrostatic  pressure 

and  will  not  be  affected  to  a  significant  amount  by  gas 

diffusion  and  temperature  effects.   A  1  mm  diameter  bubble 

at  150  m  depth  will  increase  in  size  to  2.5  mm  in  ascending 

isothermally  to  the  surface. 

A  sample  calculation  of  this  effect  is  shown  in  Table 

II  and  shows  that  this  bubble  will  shift  its  resonant 

frequency  by  several  kilohertz  during  its  ascent. 

TABLE  II 
Typical  Bubble  Size  Expansion 


RADIUS  in  mm 
(600  ft.) 


RADIUS  in  mm 
( surface) 


RESONANCE   FREQUENCY  in  kHz 
(600  ft.)    (surface) 


.200 

.5 

70.4 

16.0 

.500 

1.23 

28.2 

6.4 

1.0 

2.5 

14.1 

3.2 

2.0 

4.9 

7.0 

1.6 

10.0 

24.6 

1.4 

.13 
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E.    RISE  OF  BUBBLES 

The  rate  of  rise  of  bubbles  has  been  thoroughly 

documented.  [Refs.  22,  23,  24,  25]   The  speed   at  which  a 

bubble  rises  depends  primarily  upon  the  size  and  shape  of 

the  bubble,  and  can  be  categorized  as  in  Table  III. 

[Ref.  26:  p.  72] 

TABLE  III 

Characteristics  of  Rising  Bubbles 

SHAPE      DIAMETER  in  MM      SPEED  in  cm/s   Re  NUMBER    MOTION 

spherical         .5  5-10  <  200    rectilinear 

ellipsoidal      1-15  10-30      200<Re<4700   zig-zag 

helical 
oscillating 

spherical-        >15  30-40        >4700       erratic 

capped  unstable 

Increasing  bubble  size  causes: 

1.  increased  coefficient  of  drag 

2.  increased  vertical  speed 

3.  nonlinear  motion 

4.  oscillating  motion 

5.  wake  shedding 

Turbulence  seems  to  be  the  controlling  factor  in  the  motion 
of  bubbles  larger  than  1,200  mm.    [Ref.  27:  p.  121]   "Bubbles 
smaller  than  about  0.1  cm  in  water  rise  along  smooth  paths 
and  do  not  generate  any  measureable  sound  when  rising." 
[Ref.  28:  p.  25] 
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Surface  contaminants  affect  bubble  behavior  in  many  ways. 
These  surfactants  [Ref.  29]  will  have  their  greatest  effect 
on  the  terminal  speed  of  a  bubble  by  affecting  transition 
from  rectilinear  motion  to  oscillating  motion.   The  bubble's 
wake  structure  is  changed  causing  delayed  boundary  layer 
separation,  delayed  vortex  shedding  and  smaller  bubble  wakes. 

Terminal  speed  will  be  a  function  of  bubble  diameter, 
surfactants  and  Reynold's  number.   Approximate  increase  in 
vertical  speed  due  to  multiple  bubbles  range  up  to  a  factor 
of  1.5  times  the  single  bubble  velocities.    [Ref.  30] 

F.    BUBBLE  NOISE 

Previous  investigations  [Ref.  31]  show  that  greater  than 
90  percent  of  sound  produced  by  bubbles  in  water  occurs  during 
bubble  formation. 

12 Q.       n       o 


o 


time 


10  ms 


Figure  1.   Sound  produced  by  a  bubble  at  an  orifice 
This  sketch  [Ref.  32:  p.  23]  correlates  the  process  of  an 
individual  bubble  leaving  an  orifice  with  the  observed  time 
variation  of  the  radiant  noise.   The  peak  sound  pressure 
occurs  at  the  time  the  bubble  separates  from  the  orifice. 

Strasberg  [Ref.  33:  p. 24]  predicted  the  peak  sound 
pressure  at  distance  d  to  be. 
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iKms 


1 


K  6        ' 


where,   Q  =   total  dissipation  constant 

Q  =   .014  for  an  air  bubble  in  water 

Y  =   total  volume  rate  of  bubble  formation 

r'  =   excess  pressure 

This  works  out  to  be, 

2    132     da   re   i  ^  Fa   at     1  m. 

for  a  bubble  rate  of  10/sec  for  f  =  1.0  kHz 
After  separation,  larger  bubbles  generate  continuous 
sound  levels  of  relatively  small  amplitude  (46  to  66  dB  re 
1  ;tiPa  at  Im).  [Ref.  34:  p. 25] 

This  is  apparently  due  to  oscillations  induced  by  their 
irregular  flow  path.   Smaller  bubbles  produce  less  or 
immeasurable  sound  when  rising.   The  coalescing  of  bubbles 
cause  one  tenth  the  sound  of  bubble  formation. 

G.    PROPERTIES  OF  BUBBLE  SCREENS 

Studies  on  the  acoustic  theory  of  bubbles  and  bubble 
screens  were  done  during  World  War  II.    [Refs.  35,  36,  37] 
Efforts  in  the  United  States  were  concentrated  through  the 
National  Defense  Research  Committee  (NDRC)  under  a  broad 
research  program  supported  by  the  Office  of  Scientific 
Research  and  Development  (OSRD) .   The  prominent  NDRC 
organizations  whose  findings  were  published  by  the  U.S. 
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Navy  after  the  war  as  NDRC  Summary  Technical  Reports  were 
the  University  of  California  Division  of  War  Research 
(UCDWR) ,  the  Columbia  University  Division  of  War  Research, 
the  Columbia  University  Sonar  Analysis  Group,  and  the  Woods 
Hole  Oceanographic  Institution.   The  UCDWR  was  the  largest 
of  the  associated  groups  and  came  to  be  known  as  the  Sonar 
Data  Division,  accounting  for  a  significant  amount  of 
empirical  studies  on  ship's  wakes,  scattering  and  absorption 
of  bubbles,  and  sound  transmission  through  bubble  screens. 
Carstenson  and  Foldy  were  part  of  this  effort  and  in  1947 
they  published  their  study  of  the  propagation  of  sound 
through  a  bubble  screen.   Their  screen  was  approximately 
17  in.  long  and  10  ft.  high  with  thickness  varying  from 
3  in.  to  6  in.   The  pulsed  screen  was  analyzed  as  if 
bubbles  were  uniformly  distributed  between  two  infinite 
parallel  planes  and  the  speed  of  sound  in  the  screen  was 
assumed  to  be  approximately  the  same  as  that  of  water. 

For  the  case  of  continuous  -  flow  screens,  the  bubbles 
were  of  varying  sizes  and  speed  of  sound  in  the  layer 
considerably  different  from  that  in  water. 

Their  results  supported  the  resonant  bubble  theory  for 
pulse  screens  at  frequencies  from  15  kHz  to  35  kHz.   However 
their  data  from  5  kHz  to  45  kHz  for  continuous  flow  screens 
of  bubbles  of  varying  sizes  were  not  very  satisfactory. 
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Throughout  their  study,  the  thickness  of  the  screen  was 
assumed  to  be  the  arithmetic  average  of  the  observed  layer 
width.   The  density  of  the  screen  was  based  on  observations 
of  the  average  number  of  bubbles  present  per  unit  volume 
and  of  the  average  bubble  size. 

For  the  continuous  flow  case,  the  screen  had  a  "..fairly 
sharply  defined  core  of  larger  bubbles  ...  but  in  front  and 
back  there  is  a  gradually  tapering  distribution  of  smaller 
bubbles."   [Ref.  38:  p.SOlJ   As  a  result,  the  speed  of  sound 
was  a  gradually  changing  function  of  the  distance  normal  to 
the  layer  interface  plane  causing  reduced  reflections  of 
sound . 

Wood  and  Spitzer  [Refs.  39,40j  proposed  theoretical  sound 
speeds  which  would  occur  for  given  concentrations  of  air  by 
volume  in  water.   Figure  2  illustrates  those  results  and 

experimental  verification  by  others.   jRefs.  41,  42]   For 

-2         -3 
concentrations  between  10    and  10   ,  readily  obtainable  m 

the  laboratory  or  in  the  field,  sound  speeds  as  low  as 

100  -  300  m/sec  are  found. 

H.   APPLICATIONS  OF  LAYERS 

There  are  extensive  studies  of  the  reflection  and  trans- 
mission of  energy  in  layered  media.   [Refs.  45,  46]   Numerous 
applications  of  the  use  of  layers  range  from  reflection 
reduction  for  optical  lens  coatings  to  absorbent  materials 
for  architectural  acoustics.   The  ratio  of  layer  thickness 
to  the  wavelength  of  the  energy  in  the  layer  is  a  very  sig- 
nificant parameter  due  to  the  standing  wave  pattern  produced. 
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Figure  2.   Sound  speed  in  air-water  mixtures 
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III.    BUBBLE  LAYER  AT  CARR  INLET 

A.    OCEANOGRAPHIC  CHARACTERISTICS 

The  chart  of  Carr  Inlet  Acoustic  Range  in  Appendix  A 
indicates  a  large,  irregularly  shaped  area  surrounding  the 
test  site.   The  site  is  roughly  5  nm  long  and  3  nm  wide. 
It  is  proposed  to  consider  in  this  study  a  bubble  screen 
crossing  the  mouth  of  the  inlet,  a  distance  of  approximately 
2900  m.   This  position  places  the  screen  directly  between 
commercial  boat  traffic  to  the  southeast  and  the  test  site. 

The  water  depth  at  this  position  varies  from  approximately 
450  ft.  at  the  relatively  level  bottom  at  the  middle  of  the 
inlet,  to  relatively  steep  shore  gradients  approaching  Fox 
Island  and  McNeil  Island  as  shown  in  Figure  3. 

The  frequencies  of  interest  for  attenuation  will  be 
broadband  (100  Hz  to  10  kHz+). 

Oceanographic  studies  [Ref.  47j  at  Carr  Inlet  indicate 
that  the  tides  and  currents  which  are  likely  to  be  encountered 
will  be  variable  but  less  than  one-half  kt.   There  is  no 
significant  evidence  for  underground  springs  or  bottom 
currents  which  would  cause  any  significant  salinity  or 
temperature  gradients.   Seasonal  variations  in  water 
temperature  and  salinity  occur  due  to  the  annual  air  tem- 
perature cycles  and  rain  runoff.   However,  no  significant 
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Figure  3.   Position  of  bubble  screen  in  sensitivity  analysis. 
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impact  is  anticipated  on  the  sound  velocity  profile.   So, 
for  the  distances  encountered  in  this  application,  the 
sound  speed  in  water  will  be  assumed  to  be  constant 
(C  =  1500  m/s) . 

B.    SCREEN  POSITION 

Considering  these  oceanographic  factors,  the  bubble 
screen  will  be  assumed  to  be  produced  at  the  bottom  of 
the  inlet  and  to  rise  vertically  with  constant  width.   An 
estimation  for  actual  layer  width  spreading  will  be 
presented  in  Chapter  V. 
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IV.    SENSITIVITY  ANALYSIS 

A.    THEORY 

At  frequencies  removed  from  the  resonant  frequency  of 
the  bubbles,  a  bubble  layer  is  expected  to  reflect  sound  as 
a  result  of  two  interrelated  effects  on  the  acoustic 
properties  of  the  medium.    [Ref.  48:  P. 925]   First,  the 
compressibility  of  the  water  is  changed  with,  the  introduction 
of  the  bubbles  which  in  turn  changes  the  sound  speed  in  the 
layer.   Second,  density  will  be  changed.   For  the  bubble  con- 
centrations considered  in  this  thesis,  this  latter  change  is 
small  and  the  principal  acoustic  characteristic  of  the  bubble 
layer  will  be  its  sound  speed.   The  sound  speed  in  bubble 
layers  was  addressed  by  Wood  [Ref.  49j  wh-ose  work  was  con- 
firmed experimentally  by  Silberman.    [Ref.  50:  p.90] 

In  the  present  work,  the  acoustic  effect  of  a  bubble 
screen  will  be  analyzed  in  terras  of  the  sound  power  trans- 
mission coefficient  using  a  computer  simulation  of  the  bubble 
screen.   As  indicated  earlier,  this  computer  simulation  will 
be  based  on  Rayleigh  reflection  theory  for  an  homogeneous 
layer  of  bubbly  water  with  constant  thickness. 

Rayleigh  reflection  theory  is  based  on  a  solution  of  a 
system  of  acoustic  equations,  describing  the  conditions  that 
must  be  satisfied  at  each  boundary  of  the  layer. 
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The  boundary  conditions  at  the  interfaces  are  continuity 
of  the  acoustic  pressure  and  of  the  normal  component  of  the 
particle  velocity. 

For  oblique  incidence  and  when  the  media  in  front  and 
behind  the  layer  are  the  same,  the  sound  power  transmission 


:oefficient  C><1  reduces  to  [Ref.  .51], 

i 


^T  = 


1-  ' 


(1) 


This  coefficient  represents  the  ratio  of  the  intensity  of 
the  transmitted  sound  to  the  intensity  of  the  incident  sound 
The  geometry  of  the  problem  is  shown  in  Figure  4,  and 
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Figure  4.   Bubble  layer  reflection  geometry 
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'-^j-   is  a  function  of  width  and  acoustic  impedance  of  the 
layer  and  depends  on  the  angle  of  incidence  and  frequency  of 
the  incident  sound. 

The  subsequent  sensitivity  analysis  will  indicate  the 
theoretical  performance  expected  when  parameters  vary. 

B.    SENSITIVITY  ANALYSIS  RESULTS 

1 .   Parameter  description 

Data  were  taken  using  the  options  in  Table  IV.   A 

change  in  symbols  was  required  for  use  in  the  computer  model 

The  new  symbols  used  to  represent  the  parameters  are  also 

listed  for  the  options  in  Table  IV. 

TABLE  IV 
Analysis  Options 


OPTION 

ORDINATE 

ABSCISSA 

CONSTANT 

PARAMETER 

1 

T3 

L2 

Al 

A 

2 

T3 

Al 

L2 

A 

3 

T3 

Al 

A 

L2 

4 

T3 

L2 

A 

Al 

5 

T3 

A 

L2 

Al 

6 

T3 

A 

Al 

L2 

where , 

T3  =  sound  power  transmission  coefficient 

L2  =  layer  width  divided  by  the  wavelength  in  layer 

Al  =  incident  angle 

A    =   sound  speed  ratio  C,       /  C   ^ 

layer  /   water 
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Each  of  these  options  will  be  addressed  in  the  following 

analysis.   Values  of  A  and  their  corresponding  fraction 

volume  of  air  bubbles  are  shown  in  Table  V, 

TABLE  V 

Sound  speed  ratios  for  given  air  concentrations 

Fraction  by  volume  of 
A  air  in  water 


0.5  2  X  lO'"^ 

0.4  5  X  10~^ 

0.3  8  X  lO""^ 

0.2  2  X  lO"^ 

0.1  8  X  lO"^ 

These  values   of  A  may  reasonably  be  expected  in  a 
practical  bubble  screen. 

The  value  of  L2  was  obtained  by  incorporating  a  normal- 
ization factor  which  forced  the  data  to  repeat  itself  after 
L2  reached  a  value  of  1.0.   This  was  done  for  ease  in 
plotting  and  clarity  in  data  presentation.   The  normalization 
factor  Ll  is  found  by  noting  that  in  Eq.  (1)  , 

c<^=  1  ,       for-      i^k=  nTC  r\=  0,1,2.,... 

for  what  values  of  layer  width/wavelength  (  ■=- I  c 

/\ivi.  lauer) 


T 


I^[x  COS  a,  =  nrr 
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Let    n    =    1 , 

1    =±L        and.     bv    Snell's     Law.        COS  B,   ="Vi  "  A  S//^^ 


where    Ll    = "Y"       and,    by    Snell's    Law, 


Therefore,    Ll    =  1   (^.ijl  -  A"^  5  i/tQ^  and    L2    = 


h 


^^Ll 


2 .  The  Computer  Program 

The  WANG  SYSTEM  2200S  desk  top  computer  was  used  to 
run  the  analysis  and  provide  data  output.   The  program, 
"0PT6",  was  written  in  BASIC  language  and  provides  a  defini- 
tion of  the  variables  and  symbols  at  the  beginning  of  the 
program. 

The  plotting  option  is  selected  when  the  program  is 
executed.   Each  option  has  its  own  subroutine  which 
interrogates  the  operator  for  the  bubble  screen  parameters 
required  to  produce  the  selected  data  output.   Step  sizes 
for  all  the  parameters  are  pre-selected .   The  sound  power 
transmission  coefficient  is  plotted  vertically  against  the 
changing  option  parameters. 

3 .  Option  1 

When  the  sound  power  transmission  coefficient  is 
plotted  against  the  normalized  layer  ratio,  T3  cycles  be- 
tween maxima  and  minima.   Minimum  sound  power  is  transmitted 
through  the  layer  at  odd  integer  multiples  of  one-quarter 
wavelengths.   Figures  5  to  10  show  the  dependence  of  T3  on 
L2  of  a  fixed  Al  for  a  variety  of  values  of  A. 
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Figure  5.   Sound  power  transmission  coefficient  vs 
layer  ratio,  parameter  A  for  constant 
Al  =  0° 
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Figure  6.   Sound  power  transmission  coefficient  vs. 
layer  ratio,  parameter  A  for  constant 
Al  =  20^^ 
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Figure  7.   Sound  power,  transmission  coefficient  vs 
layer  ratio,  parameter  A  for  constant 
Al  =  40° 
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Figure  8.   Sound  power  transmission  coefficient  vs 
layer  ratio,  parameter  A  for  constant 
Al  =  60° 
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Figure  9 


Sound  power,  transmission  coefficient  vs 
layer  ratio,  parameter  A  for  constant 
Al  =  80° 
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Figure  10.   Sound  power  transmission  coefficient  vs 
layer  ratio,  parameter  A  for  constant 
Al  =  85° 
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.  .   The  following  observations  can  be  made: 

a.  Regardless  of  angle  of  incidence,  if  the  sound 
speed  ratio  is  less  than  0.1,  attenuation  of  at  least  lOdB 
can  be  expected  over  at  least  the  middle  75  percent  of  the 
layer  ratio  range. 

b.  Regardless  of  the  sound  ratio  or  the  angle  of 
incidence,  there  will  always  be  at  least  5-10  percent  of  the 
layer  ratio  which  will  have  only  negligible  attenuation. 

4 .   Option  2 

When  the  sound  power  transmission  coefficient  is 
plotted  against  the  angle  of  incidence,  the  maximum  trans- 
mission always  occurs  at  normal  incidence  and  the  minimum 

o 
transmission  at  angles  approaching  90  .   Figures  11  to  16 

are  plots  with  the  layer  ratio  held  constant  and  the  sound 

speed  ratio  varied  as  a  parameter. 

The  following  observations  can  be  made: 

a.  When  the  sound  speed  ratio  reaches  0.1,  at 
least  lOdB  attenuation  can  be  expected  over  all  angles  of 
incidence  and  over  at  least  75  percent  of  the  layer  ratio 
range. 

b.  The  only  values  of  sound  speed  which  will  give 
a  uniform  value  of  attenuation  over  all  angles  of  incidence 
are  speed  ratios  of  less  than  0.1.   This  situation  would 

require  the  bubble  screen  to  maintain  a  volume  fraction 

-2 
of  air  to  water  of  greater  than  approximately  10 
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Figure  11.  Sound  power  transmission  coefficient  vs. 
incident  angle,  parameter  A  for  constant 
L2  =  0 
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Figure  12. 


Sound  power  transmission  coefficient  vs. 
incident  angle,  parameter  A  for  constant 
L2  =  .1 
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Figure  13. 


Sound  power  transmission  coefficient  vs. 
incident  angle,  parameter  A  for  constant 
L2  =  .2 
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Figure  14.  Sound  power  transmission  coefficient  vs. 
incident  angle,  parameter  A  for  constant 
L2  =  .  3 
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Figure  15. 


Sound  power  transmission  coefficient  vs. 
incident  angle,  parameter  A  for  constant 
L2  =  .4 
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Figure  16 . 


Sound  power  transmission  coefficient  vs. 
incident  angle,  parameter  A  for  constant 
L2  =  .  5 
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5 .   Option  3 

As  in  the  plots  for  option  2,  the  sound  power 
transmission  coefficient  is  plotted  against  the  angle  of 
incidence.   This  time,  though,  the  sound  speed  ratio  is 
plotted  as  a  parameter.   See  Figures  17  to  22. 

The  maximum  sound  transmitted  occurs  at  normal 
incidence  and  the  minimum  always  occurs  as  the  angle  of 
incidence  approaches  90  . 

The  following  observations  can  be  made: 

a.  The  significance  of  a  change  in  the  sound 
speed  ratio  is  very  clear  in  this  option.   If  Figure  18 
(sound  speed  ratio  =  0.1)  is  compared  to  Figure  22  (sound 
speed  ratio  =  0.5) ,  a  dramatic  decrease  of  attenuation  is 
seen  to  occur. 

b.  At  low  angles  of  incidence,  values  expected  for 
attenuation  are  relatively  constant.   At  high  angles  of 
incidence,  attenuation  values  are  changing  rapidly, 
indicating  less  confidence  for  a  given  predicted  layer  width 

5 .   Option  4 

This  option  plots  the  sound  power  transmission 
coefficient  against  the  normalized  layer  ratio.   The  angle 
of  incidence  is  the  parameter  which  is  varied  for  constant 
values  of  sound  speed  ratios  (Figures  23  to  29) .   Again  it 
is  seen  that  values  for  the  sound  power  transmission 
coefficient  cycle  between  maxima  and  minima.   The  minimum 
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Figure  17 . 


Sound  power  transmission  coefficient  vs. 
incident  angle,  parameter  L2  for  constant 
A  =  .05 
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Figure  18 . 


Sound  power  transmission  coefficient  vs. 
incident  angle,  parameter  L2  for  constant 
A  =  .1 


48 


13  ^-0 


p^^^^^^^^^r™^ 


.5- 


A 

A1 

L2 

T3 


A1 
SOUND  SPEED  RATIO 
INCIDENT    ANGLE 
LAYER  WIDTH/ WAVELENGTH  RATIO 
OUND  POWER  TRANSMISSION  COEFFICIENT 


_     «^ 


Figure    19. 


Sound  power  transiaission  coefficient  vs. 
incident  angle,  parameter  L2  for  constant 
A  =  .  2 
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Figure  20 . 


Sound  power  transmission  coefficient  vs. 
incident  angle,  parameter  L2  for  constant 
A  =  .  3 
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Figure  21.   Sound  power  transmission  coefficient  vs. 
incident  angle,  parameter  L2  for  constant 
A  =  .4 
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Figure  22  . 


Sound  power  transmission  coefficient  vs. 
incident  angle,  parameter  L2  for  constant 
A  =  .  5 


52 


T3 


io 


A1  =  0»  TO  89» 
—     ^^-^-^  1  I  I  I  I  I  I  1 


I  I  l,i  I  I 


.5 


A 

A1 

L2 

T3 


SOUND  SPEED  RATIO 

INCIDENT    ANGLE 

LAYER  WIDTH/ WAVELENGTH   RATIO 

SOUND  POWER  TRANSMISSION  COEFFICIENT 


Figure  23  . 


Sound  power  transmission  coefficient  vs 
layer  ratio,  parameter  Al  for  constant 
A  =  .05 
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Figure  24. 


Sound  power  transmission  coefficient  vs 
layer  ratio,  parameter  Al  for  constant 
A  =  .1 
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L2 
A       =     SOUND  SPEED  RATIO 
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Figure  2  5 . 


Sound  power  transmission  coefficient  vs 
layer  ratio,  parameter  Al  for  constant 
A  =  .  2 


55 


3^-° 

V 

i 

.5- 

|i 

\    \                                                                                                           i 

/ 

\        \ 

x>-^^il=2>-</ 

• 

\ 

\      ^^^--^      ^.-^^ 

- 

\ 

\v.^                       A'b\y^ 

K^. 

8  90 

.  .  -,  -f^  . 

.5 


A 

A1 
L2 

T3 


1.0 
L2 


SOUND  SPEED  RATIO 

INCIDENT    ANGLE 

LAYER  WIDTH/ WAVELENGTH  RATIO 

SOUND  POWER  TRANSMISSION  COEFFICIENT 


Figure  25 . 


Sound  power  transmission  coefficient  vs 
layer  ratio,  parameter  Al  for  constant 
A  =  .  3 
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Figure  27. 


Sound  power  transmission  coefficient  vs 
layer  ratio,  parameter  Al  for  constant 

A  =  .4 
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Figure  28.   Sound  power  transmission  coefficient  vs 
layer  ratio,  parameter  Al  for  constant 
A  =  .5 
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Figure  29 . 


Sound  povrer  transmission  coefficient  vs 
layer  ratio,  parameter  Al  for  constant 
A  =  .9 
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values  occur  at  layer  ratios  which  correspond  to  odd 
multiples  of  one-quarter  wavelength. 

The  following  observations  can  be  made: 

a.  Changes  in  attenuation  due  to  variations  in 
the  angle  of  incidence  are  very  slight  for  angles  less 
than  approximately  45  ,  but  the  attenuation  in  uniformly 
predictable  and  by  far  the  greatest  at  angles  greater 
than  80°-85°. 

b.  Variations  in  the  sound  speed  ratio  cause  far 
greater  effects  on  the  attenuation  when  the  sound  speed 
ratio  is  greater  than  0.2. 

7  .   Option  5 

The  sound  power  transmission  coefficient  is  plotted 
against  the  sound  speed  ratio  in  this  option.   The  bubble 
layer  ratio  is  kept  constant  while  the  angle  of  incidence 
is  varied.   A  much  different  view  of  the  data  is  seen  in 
Figures  30  -  34. 

The  following  observations  can  be  made: 

a.  The  best  attenuation  occurs  at  very  low  sound 
speed  ratios  as  was  also  seen  in  options  1  and  2. 

b.  There  is  a  significant  amount  of  attenuation 
gained  for  angles  of  incidence  past  75  . 

8 .   Option  6 

In  this  option,  the  sound  power  transmission 
coefficient  is  plotted  against  the  sound  speed  ratio  with 
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the  parameter  of  layer  ratio  varied  for  constant  values 
of  angles  of  incidence  in  Figures  35  to  38. 

The  following  observation  can  be  made: 
a.   Layer  width  ratio  curves  vary  only  slightly 
indicating  that  layer  width  is  not  an  important  parameter 
In  addition,  there  is  still  approximately  10  percent  of 
the  layer  ratio  range  (L2<  .05)  virtually  unaffected  by 
sound  speed  ratio  or  angle  of  incidence. 
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A  =  SOUND  SPEED  RATIO 

A1  =  INCIDENT    ANGLE 

L2  =  LAYER  WIDTH/ WAVELENGTH  RATIO 

T3  =  SOUND  POWER  TRANSMISSION  COEFFICIENT 


Figure  30  . 


Sound  power  transmission  coefficient  vs, 
sound  speed  ratio,  parameter  Al  for  constant 
L2  =  .05 
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T3     =     SOUND  POWER  TRANSMISSION  COEFFICIENT 


Figure  31 . 


Sound  power  transmission  coefficient  vs. 
sound  speed  ratio,  parameter  Al  for  constant 
L2  =  .1 
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Figure  32. 


Sound  power  transmission  coefficient  vs. 
sound  speed  ratio,  parameter  Al  for  constant 
L2  =  .  2 
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Figure  33 . 


Sound  power  transmission  coefficient  vs. 
sound  speed  ratio,  parameter  Al  for  constant 
L2  =  .3 
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Figure  34. 


Sound  power  transmission  coefficient  vs. 
sound  speed  ratio,  parameter  Al  for  constant 
L2  =  .4 
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Figure  35.   Sound  power  transmission  coefficient  vs. 

sound  speed  ratio,  parameter  L2  for  constant 
Al  =  0° 
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A  =  SOUND  SPEED  RATIO 

A1  =  INCIDENT    ANGLE 

L2  =  LAYER  WIDTH/ WAVELENGTH  RATIO 

T3  =  SOUND  POWER  TRANSMISSION  COEFFICIENT 


Figure  36.   Sound  power  transmission  coefficient  vs. 

sound  speed  ratio,  parameter  L2  for  constant 
Al  =  3Q 
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Figure  37.   Sound  power  transmission  coefficient  vs. 

sound  speed  ratio,  parameter  L2  for  constant 
Al  =  60° 
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Figure  38.   Sound  power  transmission  coefficient  vs. 

sound  speed  ratio,  parameter  L2  for  constant 
Al  =  85 
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V.    BUBBLE  SCREEN  APPROXIMATIONS  AND  RECOMMENDATIONS 

A.  FREQUENCY  BANDS  WITH  LITTLE  ATTENUATION 

Options  1  and  6  shows  that  there  will  be  at  least  10 
percent  of  the  layer  ratio  range  virtually  unaffected  by 
a  bubble  screen.   So  for  any  given  spectrum,  there  will  be 
frequency  bands  of  small  attenuation  which  occur  when  the 
screen  thickness  is  an  integer  multiple  of  one-half 
wavelengths  in  the  layer. 

B.  LAYER  SPREADING 

Answers  to  questions  on  layer-width  spreading  were  not 
found  in  the  literature.   Forcing  parameters  which  are 
present  in  Carr  Inlet  which  would  cause  a  bubble  screen  to 
lose  its  simulated  parallel  interfaces  with  the  surrounding 
media  are  bubble  interactions,  turbulent  bubble  motion, 
and  inlet  currents.   In  order  to  understand  these  effects, 
it  might  be  better  to  conduct  an  experimental  study  at 
Carr  Inlet  and  not  a  theoretical  study  because  predictability 
rapidly  disappears  due  to  turbulent  effects  when  bubbles 
grow  larger  than  1  to  2  cm. 

However,  an  estimate  of  the  layer  spreading  may  be 
taken  from  a  linear  extrapolation  of  the  observed  laboratory 
results  of  Carstenson  and  Foldy.   [Ref.  52]   Their  screen 
spread  up  to  6  in.  for  a  rise  of  10  ft.   Figure  39  shows 
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Figure  39.   Bubble  screen  spreading  effect 
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approximate  layer  thicknesses  for  a  linear  e:xtr  apolat  ion . 
The  corresponding  screen  volume  increase  in  this  case  is  a 
factor  of  60  from  600  ft.  to  the  surface.   At  the  same  time, 
the  bubbles  formed  at  600  ft.  are  expanding  their  volume 
due  to  decreased  hydrostatic  pressure  by  a  factor  of  about 
18  atmospheres.   So  for  this  hypothetical  case,  the  fraction 
volume  of  air  to  water  changes  by  a  factor  of  about  one-third 
This  could  cause  a  change  of  speed  of  sound  of  200  -  300  m/s . 

When  considering  an  approximation  like  this,  the  number 
of  unknowns  involved  indicate  the  uncertainty  inherent  in 
projecting  any  laboratory  layer  spreading  data  to  the 
dimensions  of  Carr  Inlet. 

C.    EXPERIMENTAL  SCALING 

The  factors  to  be  considered  for  proper  scaling  of 
laboratory  experiments  are: 

1.  Pressure  change  during  bubble  ascent 

2.  Height,  length,  thickness  of  the  bubble  layer 

3.  Angle  of  incidence  and  frequency  of  incident  sound 

4.  Path  lengths  of  sound 

Of  these  factors,  the  pressure  change  is  the  most 
difficult  to  scale  because  volume  expansion  of  the  bubble 
would  cause  a  significant  change  in  the  speed  of  sound. 

On  site  experiments  at  Carr  Inlet  appear  to  be  the  only 
viable  method  to  test  the  results  of  this  sensitivity 
analysis . 
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D.    RECOMMENDATIONS  AND  ALTERNATIVES 

The  purpose  of  this  thesis  was  to  investigate  the 
benefits  of  a  bubble  screen  which  may  offer  some  isolation 
from  underwater  noise  interference  at  Carr  Inlet.   As  a 
result  of  this  sensitivity  analysis  and  further  investiga- 
tion, the  following  recommendations  are  made: 

1.   An  alternate  screen  position  as  shown  in  Figure  40 


®  TEST  SITE 

\\\  ALTERNATE  SCREEN 


Figure  40.   Alternate  bubble  screen  position 
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a.  The  advantages  of  this  position  are: 

(.1)  Only  a  fraction  of  the  piping  (one  estimate 
for  piping  to  cross  the  inlet  for  a  uniform  screen  required 
either  a  large  number  of  in  line  pressure  regulators  or 
almost  100  km  of  piping) . 

(2)  Shallower  depth  allowing  less  dispersion  of 
bubble  layer. 

(3)  Exploitation  of  the  high  attenuation  pre- 
dicted in  this  sensitivity  at  steep  angles  of  incidence. 

b.  The  disadvantages  include: 

CD  Increased  noise  interference  when  screen 
is  placed  closer  to  the  test  site. 

(2)  Low  frequency  sound  still  not  affected. 

C3)  Diffraction  of  sound  around  edges  reduces 
any  attenuation.   Further  study  for  any  perforated  pipe 
bubble  screen  mechanism  should  include  pipe  suspension  or 
anchoring  devices,  pipe  cleaning  methods  and  compressor 
requirements . 

2.  Conduct  an  on-site  layer  spreading  experiment. 

3.  As  an  alternative  to  the  bubble  screen  solution, 
investigate  the  benefits  of  adaptive  beam  forming  to  null 
out  bothersome  interference  directions. 
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VI .    CONCLUSIONS 

A.  The  insulating  capability  of  a  bubble  screen  at  Carr 
Inlet  is  a  function  of  the  angle  and  frequency  of  incident 
sound  energy  and  the  width  and  speed  of  sound  of  the  bubble 
layer.   Of  these  parameters,  the  angle  of  incidence  and  the 
speed  of  sound  in  the  layer  are  the  most  important  to  consider 

B.  There  will  be  approximately  10  percent  of  the  frequency 
spectrun  over    which  virtually  no  attenuation  of  sound  will 
occur  no  matter  what  speed  of  sound  could  reasonably  be 
acheived  in  the  bubble  layer. 

C.  The  maximum  consistent  attenuation  which  could  be 
expected  is  approximately  10  dB  at  bubble  concentrations 
greater  than  about  5  x  10 

These  conclusions  are  based  on  the  assumption  that  the 
air-water  mixture  in  the  layer  is  an  ideal,  homogeneous  fluid. 
As  a  next  step  in  the  improvement  of  this  model,  the  case 
of  a  lossy  medium  should  be  incorporated  to  include  the 
effects  of  absorption  of  a  bubble  screen. 
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APPENDIX  A:    BUBBLE  SCREEN  POSITION  AT  CARR  INLET 
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COMPUTER  PROGRAM 


10  REM  ;  PROGRAM  NAME "OPTS  " 

20  REM  ;  PURPOSE ...  PLOTS  THE  POWER  TRANSMISSION  COEFFICIENT 
(T3)  FOR  A  PLANE  WAVE  INCIDENT  ON  A  BUBBLE  SCREEN  LAYER  IN 
WATER.   THE  DESIRED  OUTPUT  DISPLAY  IS  SELECTED  FROM  ONE  OF 
THE  OPTIONS  BELOW 
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REM 
REM 
REM 
REM 
REM 
REM 
REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

INPUT 

SELECT 

IF  S  = 


DEFINITION  OF  VARIABLES 

"A"  IS  ACOUSTIC  IMPEDANCE  RATIO  =  R2/R1 

"Al"    IS  THE  INCIDENT  ANGLE  MEASURED  FROM  NORMAL 

"A2"  IS  THE  ANGLE  OF  TRANSMISSION  IN  THE  SCREEN 

"C"  IS  THE  SOUND  SPEED  RATI0=C2 / 15 00= "A"  APPROX 

"C2"  IS  THE  SOUND  SPEED  IN  THE  BUBBLE  LAYER 

"L"  IS  THE  RATIO,  BUBBLE  SCREEN  WIDTH/WAVELENGTH 

"LI"  IS  THE  NORMALIZING  FACTOR  FOR  L 

"L2"  =  L/Ll 

"T3"  IS  THE  SOUND  POWER  TRANSMISSION  COEFF (OFFICER ) 

"S"  IS  THE  SELECTED  PLOTTING  OPTION 

PLOTTING  OPTIONS 


IF 

IF 

IF 

IF 

IF 

REM 

REM 


OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
OPTION 
"OPTION 
R 

THEN 
THEN 
THEN 
THEN 
THEN 
THEN 


IS 
IS 
IS 
IS 
IS 
IS 


T3 
T3 
T3 
T3 
T3 
T3 


VS 
VS 
VS 
VS 
VS 
VS 


NUMBER  = 

300 

570 

840 

1100 

1380 

1700 


L2 

Al 

Al 

L2 

A 

A 

S 


FOR  CONST  Al  AND  PARAMETER  A 
FOR  CONST  L  AND  PARAMETER  A 
FOR  CONST  A  AND  PARAMETER  L2 
FOR  CONST  A  AND  PARAMETER  Al 
FOR  CONST  L  AND  PARAMETER  Al 
FOR  CONST  Al  AND  PARAMETER  L2 


SUBROUTINE  OPTION  1 

THIS  SUBROUTINE  PLOTS  T3  VS   L2  (Al  CONST),  PAR  A 
INPUT  "INCIDENT  ANGLE  IN  RADIANS (EA  10  DEG=.174  RAD)=",A1 
SELECT  PRINT  213(64) 
PRINT  "      INCIDENT  ANGLE=",A1 

PRINT  "      T3 (0-1.0)  VS  L2( 0-1.0),  PARAMETER  A(.l-.9)" 
SELECT  PRINT  413 
PLOT  (10)  (0,-50,"+") 
PLOT 

PLOT  (-500 
FOR  A  = 
H=0 


(10)  (50,  0,  "+") 


1  TO 


0,) 

.9  STEP  .4 
V=0:  H1=0:  V1=0 
A2  =  ARCSIN(A"SIN(A1) ) 
Ml  =  A"C0S(A1)/C0S  (A2) 
M2  =  (1/A)'-C0S(A2)/C0S(A1) 
LI  =  1/(2'^SQR(1-A'''A''-SIN(A1) 
FOR  L2  =  0  TO  1  STEP  .02 
B2  =  2"#PI"L2''-L1"C0S(A2) 


'SIN(Al))  ) 
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460  Dl  =  COS(B2)"COS(B2) 

470  D2  =  SIN(B2)''>SIN(B2) 

480  T3  =  l/(Dl+(  (M1  +  M2)"(M1  +  M2)/4'''D2)  ) 

490  H  =  INT(L2'''500)  :  V  =  INT(T3"500) 

500  PL0T((H-H1) , (V-Vl) ,HEX(FE)) 

510  HI  =  H:  VI  =  V 

520  NEXT  L2 

530  PLOT  (-H,-V, ) 

540  NEXT  A 

550  GOTO  210 

560  REM  ;  SUBROUTINE  OPTION  2 

570  REM  ;    THIS  SUBROUTINE  PLOTS  T3  VS  Al  (L  CONST),  PARAM .  A 
58  0  INPUT  "BUBBLE  LAYER/WAVELENGTH= " , L 
590  SELECT  PRINT  213(64) 

6  00  PRINT  "      BUBBLE  LAYER/WAVELENGTH= " , L 

610  PRINT  "      T3  (0-1.0)  VS  Al(0-89  DEG),  PARAM  A(.l-.9) 

620  SELECT  PRINT  413 

630  PLOT  (10)  (0,-50,"+") 

640  PLOT  (10)  (50,0,"+") 

650  PLOT  (-500,0,) 

660  FOR  A  =  .1  to  .9  STEP  .4 

670  H=0:  V=0:  H1=0:  V1=0 

580  FOR  Al  =  0  TO  1.56  STEP  .04 

690  A2  =  ARCSIN(A"SIN(A1) ) 

70  0  Ml  =  A"  C0S(A1)/C0S(A2) 

710  M2  =  (1/A)'''C0S(A2)/C0S(A1) 

7  20  B2  =  2"#PI"L"COS(A2) 
730  Dl  =  COS(B2)"COS(B2) 
740  D2  =  SIN(B2)"SIN(B2) 

750  T3  =  l/(Dl+((Ml+M2)"(Ml  +  M2)/4''-D2)  ) 
750  H=INT(A1"286)  :  V=  INT(T3'-500) 
770  PLOT  ( (H-Hl) ,(V-V1) ,HEX(FB) ) 
780  HI  =  H:  VI  =  V 

7  90  NEXT  Al 

800  PLOT  (-H,-V,) 
810  NEXT  A 
820  GOTO  210 

830  REM  ;  SUBROUTINE  OPTION  3 

840  REM  ;  THIS  SUBROUTINE  PLOTS  T3  VS  Al  (CONST  A),  PARAM.  L2 

8  50  INPUT  "IMPEDANCE  RATIO  (R2/R1)=",A 
860  SELECT  PRINT  213(64) 

870  PRINT  "      IMPEDANCE  RATIO  (R2/R1)=",A 

880  PRINT  "      T3  (0-1)  VS  Al(0-89  DEG),  PARAM  L2()-0.5)" 

890  PLOT  (10)  (0,-50,"+") 

900  PLOT  (10) (50,0,"+") 

910  PLOT  (-500,0,) 

920  FOR  L2  =  0  TO  .5  STEP  .1 

930  H=0:  V=0:  H1=0:  V1=0 

940  FOR  Al  =  0  TO  1.56  STEP  .04 

950  LI  =  1/ (2"SQR(1-A"A"SIN(A1)"SIN(A1) )) 
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960  A2  =  ARCSIN(A"SIN(A1) ) 

970  Ml  =  A"C0S(A1)/C0S(A2) 

980  M2  =  (1/A)"C0S(A2)/C0S(A1) 

990  B2  =  2"#PI"L2"L1"C0S(A2) 

1000  Dl=  COS(B2)"COS(B2) 

1010  D2  =  SIN(B2)''-SIN(B2) 

1020  T3  =  l/(Dl+(  (M1+M2)"(M1  +  M2)/U-'=D2)  ) 

1030  H  =  INT(A1"286):  V  =  INT(T3"500) 

1040  PLOT  ((H-Hl), (V-V! ) ,HEX(FB)) 

1050  HI  =  H:  VI  =  V 

1060  NEXT  Al 

1070  PLOT  (-H,-V,) 

1080  NEXT  L2 

1090  GOTO  210 

1100  REM  ;     SUBROUTINE  OPTION  4 

1110  REM  ;  THIS  SUBROUTINE  PLOTS  T3  VS  L  (CONST  A),  PARAM  Al 

1120  INPUT  "IMPEDANCE  RATIO  (R2/R10=",  A 

1130  SELECT  PRINT  213(64) 

1140  PRINT  "      IMPEDENCE  RATIO (R2/R1 )=" ,A 

1150  PRINT  "      T3(0-l)  VS  L2(0-l),  PARAMETER  Al(22  DEG  INC)' 

1160  SELECT  PRINT  413 

1170  PLOT  (10)  (0,-50,"+") 

1180  PLOT  (10)  (50,0,"+") 

1190  PLOT  (-500,0,) 

1200  FOR  Al  =  0  TO  1.56  STEP  .39 

1210  H=0:  V=0:  H1=0:  V1=0 

1220  LI  =  1/  (2"SQR(1-A"A"SIN(A1)"SIN(A1) ) ) 

1230  FOR  L2  =  0  TO  1.  STEP  .02 

1240  A2  =  ARCSIN(A"SIN(A1) ) 

1250  Ml  =  A''-C0S(A1)/C0S(A2) 

1260  M2  =  (1/A)"C0S(A2)/C0S(A1) 

1270  B2  =  2"#PI''-L2"L1"C0S(A2) 

1280  Dl  =  COS(B2)'''COS(B2) 

1290  D2  =  SIN(B2)'''SIN(B2) 

1300  T3  =  l/(Dl+((Ml  +  M2)-'(Ml  +  M2)/4''-D2)) 

1310  H  =  INT(L2"500):  V  =  INT(T3"500) 

1320  PLOT  ((H-H1),(V-V1)  ,  HEX(FB)) 

1330  HI  =  H:  VI  =  V 

1340  NEXT  L2 

1350  PLOT  (-H,-V,) 

1360  NEXT  Al 

1370  GOTO  210 

138  0  REM  ;  SUBROUTINE  OPTION  5 

1390  REM  ;  THIS  SUBROUTINE  PLOTS  T3  VS  A  (CONST  L2 ) ,  PARAM  Al 

1400  INPUT  "BUBBLE  LAYER/WAVELENGTH  RATIO=",L 

1410  SELECT  PRINT  213(64) 

1420  PRINT  "      BUBBLE  LAYER/WAVELENGTH  RATIO=",L 

1430  PRINT  "      T2(0-l)  VS  A(O-l),  PARAM  Al(22  DEG  INC) 

1440  PLOT  (10)  (0,-50,"+") 

1450  PLOT  (10)  (50,0,"+") 

1460  PLOT  (-500,0,) 
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1470  FOR  Al  =  0  TO  1.56  STEP  .39 

1480  H=0:  V=0:  H1=0:  V1=0 

1490  FOR  A  =  . 1  TO  .9  STEP  .04 

1500  A2  =  ARCSIN(A"SIN(A1)) 

1510  Ml  =  A"C0S(A1)/C0S(A2) 

1520  M2  =  (1/A)'''C0S(A2)/C0S(A1) 

15  30  B2  =  2'-#PI"L"C0S(A2) 

1540  Dl  =  C0S(B2)"C0S(B2) 

1550  D2  =  SIN(B2)"SIN(B2) 

1560  T3  =  l/(Dl+(  (Ml  +  M2)-(Ml+M2)/4''-D2)  ) 

1570  H  =  INT(A"500):  V  =  INT(T3"500) 

1580  PLOT  ((H-Hl), (V-Vl) ,HEX(FB) ) 

1590  H1=H:  V1=V 

160  0  NEXT  A 

1610  PLOT  (-H,-V,) 

1620  NEXT  Al 

1630  GOTO  210 

17  00  REM  ;  SUBROUTINE  OPTION  5 

1710  SELECT  R 

1720  REM  ;  THIS  SUBROUTINE  PLOTS  T3  VS  A  (CONST  Al) ,  PARAM  L2 

1730  INPUT"ANGLE  OF  INCIDENCE= " ,A1 

1740  SELECT  PRINT  213(64) 

1750  PRINT  "ANGLE  OF  INCIDENCE= " , Al 

1760  PRINT  "T3(0-l)  VS  A(O-l),  PARAMETER  L2(0-.5,  .1  STEPS)" 

1770  PLOT  (10)  (0,-50,"+") 

1780  PLOT  (10)  (50,0,"+") 

1790  PLOT  (-500,0,) 

1800  FOR  L2  =  0  TO  .5  STEP  .1 

1810  H=0:  H1=0:  V=0:  V1=0 

1820  FOR  A  =  .1  TO  .9  STEP  .04 

18  30  LI  =  1/(2"SQR(1-A"A"SIN(A1)"SIN(A1) ) ) 
1840  A2  =  ARCSIN(A"SIN(A1) ) 

1850  Ml  =  A"C0S(A1)/C0S(A2) 

1860  M2  =  (1/A)"C0S(A2)/C0S(A1) 

187  0  L  =  L2-L1 

1880  B2  =  2"#PI--L"C0S(A2) 

1890  Dl  =  COS(B2)''-COS(B2) 

1900  D2  =  SIN(B2)"SIN(B2) 

1910  T3  =  l/(Dl+((Ml+M2)'MMl+M2)/4''-D2)  ) 

1920  H  =  INT(A"500):  V  =  INT(T3"500) 

1930  PLOT  ((H-Hl) , (V-Vl) ,HEX(FB)) 

1940  HI  =  H:  VI  =  V 

1950  NEXT  A 

1960  PLOT  (-H,-V,) 

1970  NEXT  L2 

1980  GOTO  210 

1990  END 
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